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■  ABSTRACT  0, 

Derivations  and  program*  are  presented  for  evaluation  of  t^he 
required  input  signal -to-noi:  e  ratio  for  specified  output  performance 
of  a  Greatest-Of-device  foil  (wed  by  an  integrator,  as  a  function  of : 

N,  the  number  of  channels  in  the  Greatest-Of;  M,  the  number  of  independent 
samples  in  the  integrator;  PJ,  the  false  alarm  probability;  and  Pw,  the 
detection  probability.  Several  numerical  examples  are  presented,  for 
the  special  case  of  P^  =  0.5.  ^ 

svb  P 
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INTRODUCTION 

In  an  effort  to  minimize  equipment,  searching,  or  display  require¬ 
ments,  a  number,  N,  of  channels  of  data  are  occasionally  represented, 
instead,  by  their  largest  member.  Hopefully  in  this  fashion,  detection 
of  signal  presence  in  a  particular  channel  will  still  be  accomplished 
reliably  because  the  largest  member  will  exceed  a  threshold.  However, 
for  N>1,  the  threshold  must  be  raised  in  order  to  maintain  an  acceptable 
false  alarm  probability;  this  effect  makes  signal  detection  less  likely. 
Also,  for  low  signal-to-noise  ratio  in  one  channel  at  the  input  to  the 
Greatest-Of  device,  virtually  a  random  selection  takes  place  in  the 
Greatest-Of  device;  whether  the  integration  that  follows  can  make  up 
for  this  non-linear  action  is  not  clear. 

Here  we  will  derive  the  equations  that  characterize  performance  of 
the  Greatest-Of  device  followed  by  an  integrator,  in  terms  of  detection 
and  false  alarm  probabilities.  We  also  will  present  BASIC  programs  for 
the  HP9830A  calculator  so  that  the  reader  can  investigate  numerical 
examples  of  his  own  interest.  Some  related  past  work  is  cited  in 
reference  1. 


DERIVATIONS 


«  T 


(jreufcwt-Of 


Figure  1.  Greatest-Of  Device  Followed  by  Integrator 

The  system  of  interest  is  depicted  in  figure  1.  Input  channels 
Xj,...xN  are  statistically  independent  and  contain  either  (a)  no  signal 
in  any  channel,  or  (b)  signal  in  one  channel  only.  The  probability 
densities  of  x  under  these  two  hypotheses  are  demoted  by  p#(x)  and 


i 
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p,  (x)  respectively.  The  integrator  in  figure  1  sums  over  M  independent 
samples  and  subjects  the  output  to  a  comparison  with  a  threshold  T. 

The  Greatest-Of  device  output  is 

J  *  HHftxfa,,...,**).  (I) 

Then  the  cumulative  distribution  of  y  for  no  signal  present  is 

Q*(Y)  *  TV*V(^<Y)*».  Signal) 

=  Pk»V(xI|..i*„<  Y|*  (2) 

where  "ft (  )  is  the  cumulative  distribution  corresponding  to  density 
ym  (  ).  Therefore,  the  density  of  y  for  no  signal  present  is 

^  = N  •  (3) 

In  a  similar  fashion,  the  distribution  of  y  for  a  signal  present 

<?,(y')=  (y<  y| 

=  Fr»t  (x.,-,  \ «  y|  s;3«i)  =  Or)  (+) 

since  only  one  channel  contains  a  signal.  The  density  of  y  is 

WOT  (s)+Tf’(a)p>b>.  & 

In  order  to  evaluate  the  statistics  of  the  integrator  output  z,  we 
assume  that  M»  1,  thereby  making  z  approximately  Gaussian.  Since 


we  have 

-MJS  ,  O) 
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Therefore  the  densities  of  z,  for  no  signal  and  signal  respectively, 
are 


where 

<t>w  *  (?t)  ft) 


where 

$(*)  s  (2t)  ex  p£-  x/i),  A) 

and  where  we  define  for  future  use 

3>(x)  =  (  M 

-«o 

The  probability  of  detection  of  the  system  in  figure  1  is  given  by 

using  (8)  and  (10).  The  false  alarm  probability  is 

■  w 

We  now  specialize  to  the  case  where  we  desire 

*  0.  s.  fa) 

In  this  case,  (11)  dictates  that  the  threshold  be  set  at 

Ts^  ^  (H) 

Then  the  false  alarm  probability  (12)  becomes 

k'W*  <W 
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For  a  specified  performance  (i.e.,  specified  Pf  and  Pn) , 
(15)  can  be  solved  for  r 


=  -  £  (PP)  ir  »as. 


The  left  side  of  (16)  is  the  (linear)  deflection  criterion  at  the 
output  of  figure  1. 

In  order  to  complete  the  analysis,  we  must  now  relate 
and  8^  to  the  statistics  of  y  given  in  (3)  and  (5).  we  have 

=  ^ O'7) 

and 

j  jj [(*“>)?  lVr.lMbht' (*) .  (,ff) 


At  this  point,  we  make  the  assumption  that  the  density  functions  p 
and  p.  of  the  inputs  are  Gaussian,  as  they  might  be,  for  example,  Tf 
the  variables  x.,...,  xN  themselves  were  sums  of  a  number  of  independent 
samples.  That  Ts, 


which  are  depicted  in  figure  2. 


The  standard  deviations 
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Figure  2.  Density  Functions  of  Input 


on  both  hypotheses  have  been  assumed  equal;  this  could  be  easily 
generalized  if  desired. 

Substitution  of  (19)  in  (17)  yields  (after  some 


labor) 


/»».  *  ”•>  +<r'a* 


where 

x  W  f  .  (2') 


The  integral  in  (21)  is  a  function  only  of  N  and  need  be  numerically 
evaluated  only  once*;  as  special  cases  we  have 

«,">>  a.' £  >«>  =  #-•  t2i) 

Substitution  of  (19)  in  (18)  yields 

/Si*  "o  +  <TC*(r), 


*  A  BASIC  program  for  the  HP9830A  calculator  for  the  evaluation  of 
is  presented  in  appendix  A. 
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CN(r)  =(N-i)f^x x  +!j A*  x 

and  where  we  have  defined 


W,-  w>» 


which  is  a  (linear)  deflection  criterion  at  the  system  input, 
special  cases  of  (24),  we  have 

C,(r)  =  r?  C*(r)  =  e^V r£  (jfr)  }  Cw(o)  *  QN 


Combining  (7),  (20),  and  (23),  there  follows 

Finally,  since  from  (7) 


we  must  evaluate 


Sz=  i  I*  si3~il  'K  ■ 

But  using  (3),  (19),  and  (10),  we  find 

’’'a-1} =  J  d* 

1-r*  bH 
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(2*) 

w 

As 

•  W 

(2D 

(2») 

(21) 

(*) 
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This  equation  gives  an  explicit  expression  for  the  minimum  deflection 
r  (defined  in  (25) )  required  for  specified  performance  fF>R>*.5, 
as  a  function  of  N  and  M.  It  requires  the  numerical  evaluation  of 
°i»,  S*(r)  as  given  in  (21),  (31),  and  (24)  respectively. 

The  left  side  of  (35)  can  be  re-arranged  into  the  form 

* S x f*'1  fa) f  M £(x)$> fr-r)- N <f&) (39 

and  has  been  adopted  for  use  here. 

A  short  table  for  the  constants  and  ^  is  given  below. 


N 

a 

L 

1 

0 

1 

1 

2 

0.564190 

1 

.825645 

3 

0.846284 

1.275664 

.747975 

4 

1.029375 

1.551329 

.701224 

5 

1.162964 

1.800020 

.668980 

6 

1.267206 

2.021739 

.644924 

7 

1.352178 

2.220304 

.626033 

8 

1.423600 

2.399535 

1 

.610653 

Table  1.  Constants  and  in  (21)  and  (31). 

EXAMPLES 

For  the  example  of* 


15000,  f?  r  id *  0*4  I0*4,  (31) 

a  tabulation  of  the  right  side  of  (35)  follows.  The  numerical 

*The  inverse  -function  of  (10)  is  tabulated  in  reference  10,  pages 

976-7. 
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Table  3.  Required  Values  of  r  for  M  =  15000,  T, 
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A  plot  of  the  required  SNR  is  given  in  figure  3,  where  the 
conversion  to  dB  is  via 


i  , 

20  |o|  r  =  2o  Icy  F  > 


since  the  inputs  to  the  Greatest-Of  device  are  voltages.  It  is 
observed,  for  example,  that  13  dB  additional  signal-to-noise  ratio 
is  required  at  the  Greatest-Of  input  for  8  channels  rather  than  1. 
The  dependence  on  N  in  figure  3  is  well  approximated  by  14.3  log  N. 
There  is  a  3.4  to  3.75  dB  difference  in  the  two  curves  in  figure  3. 


To  ascertain  how  the  curve  in  figure  3  depends  on  M,  a  plot  for 
M  =  1000  and  15000  is  presented  in  figure  4  for  Pp  »  I0~3  .  The 
required  values  of  r  are  given  in  table  4. 


N 

Retired  r  hr  Pp  =  |0  * 

I 

.  CHll 

i 

.  15-46 

3 

.ion 

4 

,1426 

S 

,1719 

<o 

.3116 

7 

,>m 

9 

.369? 

TftMe.  4*.  'Re^unr*!  Vol«fj  r  -jor  Mm  ,0O0/Tp  *  0. S' 


14 


TM  No. 
TC-13-75 


The  plot  in  figure  4  for  M  =  1000  is  well  approximated  by  the  dependence 
13  log  N.  The  difference  in  the  two  curves  is  10  log  ."'=  11.8  dB  at 

N  =  1,  and  decreases  to  10.2  dB  at  N  =8. 

A  comparison  of  the  required  signal-to-noise  ratio  at  N=8  versus 
that  at  N=1  is  presented  in  table  5  and  figure  5  as  a  function  of  M 

from*  M=10  to  15000.  Thus,  less  than  5  dB  difference  exists  at  M  =  10. 


M 

it  4« 

of  (15)  4rN«l 

*iJe 

5)  WH** 

■Required  r 
■for  Ns  * 

10 

.Vnz 

.  5167 

1.704 

do 

.5642 

.3445 

1.240 

lOo 

.3010 

.  1**7 

,?5»3 

3oo 

.n«4 

.1010 

5*  31 

/o  00 

,01172 

. 05467 

.3  689 

4ooo 

.04**6 

.021*4 

.  206+ 

15-000 

.  02523 

.0  1  54  1 

.  II  37 

Table  5.  Required  Signal-to-Noise  Ratio  for  Pp  =  .5,r^*f0 


*  The  Gaussian  approximation  for  z  in  figure  1  is  being  stretched  at 
the  left  end  of  figure  5. 
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APPENDIX  A 

Program  for  (21)  and  (31) 


The  program  below  computes  aN  in  (21).  In  order  to  compute  C# 
in  (31),  line  230  is  changed  by  replacing  the  first  X  on  the  right 
side  by  Xt 2  .  Line  10  is  for  inputting  N;  lines  20-220  is  Simpson's 
rule  for  integration  with  automatic  halving;  line  230  is  the  integrand 
of  (21);  and  lines  240-340  are  for J  of  (10). 

18  N<: 


•i  Si--'  } 


28 

j  t  —  N  97 

R  ■:  3  *  O'  I ) 

3  9 

**■»  \  -■■■  -■  3 

L~, 

'•  r;i  !, 

f !  !  !:  i:  TIL  OF  300  LIVE 

6  R 

‘“’R  I  j  )  | 

?h 

PRINT 

H  =  '  H 

•  O 

33= < FRO 

l  0  >  +  F  H  S  <  3  2  > )  *  0 , 5 

90 

34  =  0 

1  1  0 

Sh=;cv 

-S  1  1*0  F, 

120 

37  =  2  .,-2 

1  30 

FOR  83 

=1  TO  05- :  STEP  2 

140 

s  4  :t-  3  i- 

F  NO  ( 0 1  +86*08 ') 

150 

NEXT  9 

3 

160 

NR  I  T i 

1  J 1  7 1 j v  83+2*104  >  *£ 

!  70 

FOR Hu  1 

£20.  1  1  j  F  ].  0 .  0 

180 

S3=S‘9  ■* 

34 

1  90 

34=8 

2  0  0 

c  IT ...  I.:" . 

? 1  0 

3  r,  •  ~  3  t  : 

0  ,  F 

2  2  0 

60  TO  1 

3  n 

2  3  O 

DEF  FN 

0  l  3  ;■  -C  1  '  ;  ;*ENr(  -0  „  F* 

•M  0 

'LIEF  r  N 

r  ( po  .> 

278 

P  R  B 

P  0 

2  6  9 

IF  P  :  7 

THEN  298 

278 

P  =  0 

288 

0  0  T  0  3 

2  0 

290 

p= i / ■ 1 

+0 > 23 1 64 i9*P  > 

308 

P---P  + 0 

.  3 1938 1 53-P* < 0 . 3565 

3 1  9 

P=P*EX 

P  -  0.5*  P  0 1 2  /  0  0  P  ?  * 

3  2  0 

IF  F'0< 

8  THEN  340 

3  3  0 

P  =  1  —  F' 

3  4  0 

RETURN 

S  r  i  i»j 

END 

I  .  ■  C.-  tr 


82-P*< 1 . 78 1477537 . p* 

2  j  5  9  ?  8  -  P  *  1 , 3 3  0 2 7  4  4  3 )  >  )) 


S, 

.5 

h 


'h. 
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APPENDIX  B. 

Program  for  (36)  — 

The  program  below  computes  r)  in  (36).  Lines  10-20  are  for 
inputting  r  and  N;  lines  30-230  is  Simpson's  rule;  lines  240-330  is 
the  integrand  of  (36);  and  lines  340-440  are  for  $  of  (10).  )iw*j 
240-340  t «  A  •£>**  ?. 


10  F 

-0. 1 137 

20  H 

::::  @ 

30  C 

1  =  1  /  S  G  R  (  2 

,,P  j  > 

40  £ 

1  =  --  8 

50  :::: 

p  --  C; 

60  F 

EM  .INTEGRAL  OF  SCO  OVER  i 8 1 i  S ?  1 

70  PRINT 

80  F 

PINT  "N  = 

"  N  j  "  R  = "  R 

90  £ 

3  ~  ( F  H  S  ( 3 1 

)  +  F  N  3  ( 8  £  >  >  *  0 . 5 

1  00 

S  4  =  0 

110 

9  5  ~  2 

120 

3  6  =  ( S  2  S 1 

)  *  0 . 5 

130 

C  "7  —  •“>  /  •” 

140 

FOR  88=1 

TO  85-1  STEP  2 

150 

S4=S4+FHS 

3 1  +  b  6  *  S  8  .) 

1  6  0 

NEXT  S3 

1  70 

MR  HE  (15 

j  1  tl  0  .1'  (  S  b  +  £  *  8  4  )  *  b  6  #  S  7  i  8  5 

189 

1  Qfi 

FORM  FIT  E2 

0  .  1  1  J  F  1  0 . 0 

200 

9  4  ~  0 " ' 

210 

r\  Pi  “  9  pi  p 

220 

86=86*0 . 5 

230 

GOTO  140 

240 

DEF  FHS  ( >■ 

> 

£50 

T1=FNP<X) 

£  b  0 

T  2  =  E  X  P  ( -  0 

.  5  *  X  t  £ ) 

270 

T  3  =  X . R 

280 

IF  N=2  THEN  310 

296  T 4  =  T 1 T <  H-£> 

360  GOTO  320 
310  T4= 1 

3 2 0  8 = C 1  * X * T4*<  <H-1 )  * T2*F N P < T 3 >  +  T1*E X P < -0.5 * T 3 1 3 - 1-| * j p T  I 
330  RETURN  S  "  . 
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